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^' Abstract 



We simulate SU{2) gauge theory with six massless fundamental Dirac fermions. We mea- 
sure the running of the coupling and the mass in the Schrodinger Functional scheme. We 
observe very slow running of the coupling constant. We measure the mass anomalous dimen- 
r^ • sion 7, and find it is between 0.135 and 1.03 in the range of couplings consistent with the 

existence of an IR fixed point. 
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1. Introduction 



Phenomenologically viable models of technicolor can be built that are based on the exis- 
tence of gauge theories with an infrared fixed point (IRFP) [ll,|2|. The latter are asymptotically 
free theories where the dynamical effects of the fermion determinant induce a non-trivial zero 
of the beta functions at low energies, leading to scale-invariance at large distances. In partic- 
ular the existence of a large mass anomalous dimension at the fixed point has been advocated 
as an important ingredient for model-building. 

The phenomenology of strongly-interacting electroweak symmetry breaking has been sum- 
marized in several interesting reviews (see e.g. Ref. |3| for early results, and Refs. [4, |5| for 
the more recent developments.). 

Several candidate theories have been singled out by analytical studies based on approxi- 
mations of the full nonperturbative dynamics. A putative phase diagram which summarizes 
nicely the most appealing options was discussed in Ref. p], see Ref. [7| for a review of recent 
results. These seminal results have triggered a number of numerical studies. Numerical sim- 
ulations of theories defined on a spacetime lattice are indeed a privileged tool to study the 
nonperturbative dynamics of these theories from first principles. Current studies have focused 
on SU(2) with two adjoint Dirac fermions |8l-ll8|. SU(3) with 8,102 fermions in the funda- 



19|-|32|. SU(3) with two sextet fermions 28|, l33l44l|. using a variety 



mental representation 

of methods. These studies have already revealed an interesting pattern of results about the 



phase diagram of strongly-interacting gauge theories, which can provide useful comparisons 
with (and hence guidance for) the analytical results. 

In this work we focus on an SU(2) gauge theory with six flavors of Dirac fermions in 
the fundamental representation, which is supposed to be close to the lower boundary of the 
conformal window for the SU(2) color group. For example, the ladder approximation predicts 
that the conformal window begins at Nf = 7|| [61], the all-orders beta function conjecture 
predicts the lower boundary of the conformal window is at Nf = 22/(2 + 7*), which depends 



on the fixed point value of < 7* < 2 |42|, |43| , and metric confinement predicts the conformal 



window begins at Nf = 6.5 [4^. This model is particularly relevant as it may be used for a 



practical realization of conformal technicolor theories |45|, |46 . 



48 



we compute 



Using the Schrodinger Functional (SF) formulation of the theory J47|, 
the running coupling and fermion mass as a function of the energy scale in the SF scheme, 
thereby deriving information on the beta function and the mass anomalous dimension. We 
find clear evidence that the running of the coupling is rather slow, and indeed compatible 



with the existence of a fixed point. As pointed out in previous studies 18[, the identification 



of a fixed point by numerical techniques is intrinsically difficult: in the vicinity of the fixed 
point the coupling changes very slowly as a function of the scale; in order to detect a slow 
running, and to be able to identify precisely the location of the fixed point, great care must 
be exercised in taming the systematic errors that arise in numerical simulations. In particular 
we need to assess critically the systematic errors that are involved in our actual procedure, 
and their propagation in the analysis of the lattice data. 

Contrary to the case of the gauge coupling, the running of the fermion mass does not slow 
down at the fixed point, and can be more easily identified by numerical methods. Results for 



the step-scaling function for the scale dependence of the renormalized mass, ap, yield a bound 
on the mass anomalous dimension at the fixed point, 7*. Currently the main source of error 
in the determination of 7=,, comes from the uncertainty in the value of the gauge coupling at 



the fixed point. First results for the anomalous dimension were obtained in Ref. [18[ for the 
SU(2) gauge theory with adjoint fermions. Rece ntly new results have been obtained for the 
SU(3) gauge theory with sextet fermions in Ref. [4l|. 

The method used and the observables considered in this work are the same as the ones 



we implemented in Ref. [18|. They are briefly summarized for completeness in Section [21 
together with the parameters of the runs that have been used for this analysis. The running 
of the coupling is encoded in the step scaling function a{u); our results for the latter are 
presented and critically discussed in Section [31 Finally the running of the mass is studied 
in Section [H the data for the mass step scaling function ap{u) compare favourably with the 
one- loop perturbative prediction, a feature that we also observed in our study of the SU(2) 
gauge theory with adjoint fermions. Even though we are unable determine whether a fixed 
point is present, our data are sufficiently precise to yield an upper bound on the value of the 
anomalous dimension throughout the range of couplings that we measure. 

2. SF formulation 



We measure the running coupling using the Schrodinger Functional method 47|, |48| . We 



follow the same procedure as in Ref. ISj except for the change from two flavours of adjoint 



fermions to six flavours of fundamental fermions. Here we briefly describe the method; for a 



full description see Ref. [18 



The Schrodinger Functional coupling is deflned on a hypercubic lattice of size L. The 
boundary conditions are chosen to impose a constant background chromoelectric fleld, and 
depend on a parameter 77. The coupling constant is then deflned as 

with k = — 24L^/a^sin(a^/L^(7r— 2?])) chosen such that 'g^ = Qq to leading order in perturbation 
theory. This gives a non-perturbative deflnition of the coupling which depends only on L and 
the lattice spacing a. We then remove the lattice spacing dependence by taking the continuum 
limit. 

We determine the mass anomalous dimension 7 from the scale dependence of the pseu- 
doscalar density renormalisation constant Zp. This is deflned as a ratio of correlation func- 
tions: 

ZpiL) = v^//p(L/2) , (2) 



as in Ref. 18 



We use the Wilson plaquette gauge action, together with fundamental Wilson fermions, 
and an RHMC algorithm with 4 pseudofermions. 

We run at Kc, deflned as the value of k for which the PCAC mass am vanishes. We measure 
am for 5 values of k for each /3 on L = 6, 8, 10, 12 lattices and interpolate to flnd k^ for each. 
We then extrapolate in a/L to determine Kc for the L = 14, 16 lattices. 



p 


L=6 


L=8 


L = 10 


L=12 


i== 14 


L^16 


2.0 


0.151788 


0.150970 


0.150576 


0.150491 


0.150334 


0.150252 


2.2 


0.147447 


0.146939 


0.146755 


0.146782 


0.146615 


0.146565 


2.5 


0.143209 


0.142825 


0.142767 


0.142811 


0.142730 


0.142716 


3.0 


0.138869 


0.138684 


0.138651 


0.138562 


0.138523 


0.138493 


3.5 


0.136130 


0.136143 


0.136104 


0.136103 


0.136096 


0.136091 


4.0 


0.134394 


0.134350 


0.134353 


0.134339 


0.134332 


0.134327 


5.0 


- 


0.132142 


0.132142 


0.132142 


0.132142 


0.132142 


6.0 


0.130753 


0.130737 


0.130748 


0.130740 


0.130739 


0.130738 


8.0 


0.129131 


0.129145 


0.129167 


0.129172 


0.129177 


0.129182 



Table 1: Values of /3, L, n used for the determination of g^ and Zp. The entries in the table are the values of 
Kc used for each combination of /3 and L. 

In practice we achieve \am\ < 0.005. At some values of /3 and L we have additional results 
at moderately small masses of \ain\ ~ 0.01, and we observe no mass-dependence within our 
statistical errors, confirming that any residual finite-mass errors are extremely small. 

2.1. Lattice parameters 

We have performed two sets of simulations in order to determine 'g^ and Zp. We use more 
values of L (six instead of four) compared to our previous simulations to improve the quality 
of the continuum limit extrapolations, and increase the step scaling factor from s = 4/3 to 
s = 3/2 to improve the measurement of the slow running of the coupling. 

To ensure our results are not affected by the presence of a bulk transition, we have measured 
the average plaquette for a range of values of (3 and k on 6'^ lattices with SF boundary 
conditions. There is a clear jump in the plaquette at low (3, suggesting the presence of a 
bulk transition. However, this disappears around (3 = 1.6. Since the lowest (3 we use for our 
measurements of 'g^ and Zp is /3 = 2.0, our results should not be affected by this transition. 

The parameters of the runs are summarised in Table [1] The values of Kc are obtained from 
the PCAC relation as described above. 

3. Results for the coupling 

We have measured the coupling in the Schrodinger Functional scheme, ^^(/3, L), for a range 
of /3, L. Our results are shown in Tabled We see immediately that the coupling is very similar 
for different L/a at constant /3, so it runs slowly. 

To analyse the running of the coupling we first define the lattice step-scaling function. 



E(m, s, a/L) = /(^o, sL/a) y^^g^j^/^^^^ 



(3) 



and its continuum limit: 



a{u,s)= lim J]{u,s,a/L) 

a/L-i-O 



(4) 



where in both cases we will use only s = 3/2. We calculate T.{u,s,a/L) from our data as 
follows: 



/S 


L=6 


L=8 


L = 10 


L^12 


L== 14 


L^IQ 


2.0 


4.941(61) 


5.521(143) 


6.053(418) 


6.109(289) 


5.913(362) 


5.726(485) 


2.2 


3.755(32) 


4.025(70) 


4.390(158) 


4.506(345) 


4.279(233) 


4.379(252) 


2.5 


2.973(21) 


3.038(37) 


3.103(72) 


3.170(67) 


3.187(174) 


3.316(151) 


3.0 


2.123(10) 


2.173(20) 


2.150(37) 


2.291(90) 


2.336(55) 


2.338(75) 


3.5 


1.660(8) 


1.707(37) 


1.730(20) 


1.751(29) 


1.825(50) 


1.715(46) 


4.0 


1.376(4) 


1.390(8) 


1.425(16) 


1.399(30) 


1.420(19) 


1.445(31) 


5.0 


- 


1.033(3) 


1.054(7) 


1.050(9) 


1.063(15) 


1.041(16) 


6.0 


0.814(1) 


0.822(3) 


0.823(6) 


0.842(6) 


0.829(12) 


0.827(11) 


8.0 


0.576(1) 


0.581(1) 


0.575(3) 


0.586(3) 


0.585(6) 


0.593(6) 



Table 2: The entries in the table are the measured vahies of g for each combination of (3 and L. 

We first discard the L = 6 data since we have found it has large lattice artifacts. We 
then interpolate the remaining data quadratically in a/ L at each /3 to find g^f/J, L) at L = 



9|, 10|, 15. Then for each L we interpolate in f3 using the functional form [18|, |32| 



1 _ (3 

t{f3,L/a) 2N 



^ f2N\' 



(5) 



We choose the smallest n which results in a x^ such that the fit is not ruled out at a 95% CL, 
and also use ra + 1 as the next best fit; this gives a 2-5 parameter fit in each case. The number 
of parameters we use for each L/a and the x^/dof for each fit are shown in Tables [3] and HI 

We now calculate S(m, s, a/L) using the fits from Eq. [5]for L = 8, 9|, 10, 10| and s = 3/2. 
Finally we extrapolate to a/L = to obtain a{u, s). Note that this extrapolation only makes 
sense when we are on the weak-coupling side of any IRFP. 

We carry out a constant continuum extrapolation, using the data at the two values of 
a/L closest to the continuum limit. We estimate the errors using a multistage bootstrapping 
procedure. The results for a{u) using the constant continuum extrapolation are plotted in 
Figure [11 where the statistical errors only are in black and the error from changing the number 
of fitting parameters are in grey. Our results are consistent with a fixed point in the region 
'g^ > 4.02. They are also compatible with the possibility that there is no fixed point at all in 
the range of couplings we have measured. However, it is clear that cr(M) is considerably below 
the 1-loop prediction at strong coupling. 

We have also carried out a linear continuum extrapolation, using the data at the four 
values of a/L closest to the continuum limit, but it is not as well constrained by our data. 
However, it is also below the 1-loop prediction at strong coupling, and has a{u)/u < 1 in the 
region 'g'^ > 4.08. This is consistent with the presence of a fixed point in the region g"^ < 4.08. 

There is a large difference between the constant and continuum extrapolations, showing 
that the systematic error from the choice of continuum extrapolation is large. This means we 
cannot determine conclusively whether or not there is a fixed point in the range of couplings 
covered by our data. 

In the vicinity of a fixed point at a coupling g*, the beta function is linear in the coupling. 



m = (3*{9-9n + - 



(6) 



t 


L/a 




8 




10 


loi 


12 


14 


15 


16 


params 
xVdof 


3 
1.92 


4 
0.54 


4 
1.24 


4 
0.48 


3 
1.66 


2 
1.54 


2 

1.88 


2 
1.36 



Table 3: Interpolation best fit parameters for g 



tt2 



t 


L/a 




8 


^3 


10 


10| 


12 


14 


15 


16 


params 
xVdof 


4 
1.25 


5 
0.58 


5 
1.42 


5 
0.54 


4 
1.19 


3 

1.61 


3 

1.06 


3 

1.05 



Table 4: Interpolation next-best fit parameters for g^ 
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Figure 1: a{u) using a constant continuum extrapolation of the two points closest to the continuum. Statistical 
error using the optimal fit parameters in black, systematic error from using different numbers of parameters 
in the fits in grey. 



/3 


L=6 


L=8 


L== 10 


L=12 


L== 14 


L=16 


2.00 


0.26636(249) 


0.27219(306) 


0.27117(241) 


0.25956(527) 


0.24564(414) 


0.24130(578) 


2.20 


0.33220(167) 


0.32060(216) 


0.30788(537) 


0.30929(137) 


0.29792(246) 


0.29198(215) 


2.50 


0.37504(32) 


0.36203(49) 


0.35095(87) 


0.34672(73) 


0.34118(88) 


0.33255(162) 


3.00 


0.40488(21) 


0.39186(31) 


0.38451(52) 


0.37955(50) 


0.37453(53) 


0.37170(56) 


3.50 


0.42102(30) 


0.40981(82) 


0.40383(32) 


0.39832(43) 


0.39461(62) 


0.39241(93) 


4.00 


0.43105(14) 


0.42192(21) 


0.41691(31) 


0.41256(34) 


0.40997(29) 


0.40746(36) 


6.00 


0.45368(8) 


0.44908(12) 


0.44597(16) 


0.44417(10) 


0.44232(15) 


0.44045(20) 


8.00 


0.46540(5) 


0.46229(7) 


0.46005(10) 


0.45822(7) 


0.45683(10) 


0.45575(9) 



Table 5: The entries in the table are the measured vahies of Zp for each combination of /3 and L. 



where /3* is a scheme-independent coefficient, which, as described in Ref. [49|], yields further 
information on the physics of these theories. In terms of the step-scahng function o"(m, s), this 
gives: 

^{u,s)=g*^{^-g*)s-^\ (7) 

We have estimated /9* by fitting cr(M, s) in the vicinity of the fixed point, and find /3* = 
0.62(12)]'^285 where the first error is statistical and the second is systematic, for those fits 
where we see a fixed point in the range of couplings covered by our data. This does not 
include the systematic error due to the choice of a constant rather than a linear continuum 
extrapolation. Better data would be needed to make the systematic errors on /3* more robust. 



4. Running mass 

We have measured the pseudoscalar density renormalisation constant Zp for a range of 
/3, L. Our results are shown in Tableland Figure [2l We see that Zp decreases with increasing 
Lja at constant /3, indicating a positive anomalous mass dimension, but the running appears 
to be slow. 

To extract 7 we ffist define the lattice step-scaling function for the mass. 



Ep('u, s, ajV) 



Zp{gQ,sL/a) 



Zp{go,L/a) 



g'^{go,L/a)=u 



and its continuum limit 



ap[u,s 



lim Sp(m, s, a/L) . 

a— >-0 



(9) 



Again, we use only s = 3/2. To calculate Tip{u,s,a/L), we proceed similarly as for 
T^{u,s,a/L). We first discard the L = 6 data, and then interpolate quadratically in a/L to 
find Zp(/3,L) at L = 9|, 10|, 15. Then for each L we interpolate in /3 using the functional 
form llSjl 



Zp(/3,L/a) = |^QQ 



(10) 



We choose the smallest n which results in an acceptable x^y ^is for the g^ fits; this gives a 5-6 
parameter fit in each case. We also use n + 1 as a next-best fit to estimate the systematic 
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Figure 2: Data for the renornialisation constant Zp as computed from lattice simulations of the Schrodinger 
functional. Numerical simulations are performed at several values of the bare coupling f3, and for several 
lattice resolutions L/a. 
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Figure 3: ap(u) using both a constant continuum extrapolation of the two points closest to the continuum, and 
a linear continuum extrapolation. Statistical error using the optimal fit parameters with a linear continuum 
extrapolation in black, systematic error including the choice of continuum extrapolation in grey. 



errors from the choice of n. The number of parameters we use for each L/a and the x^/dof 
for each fit are shown in Tables E] and [71 

We can now calculate T,p{u,s,a/L) using Eq. [H] and the fits from Eq. [101 and finally 
extrapolate to the continuum limit to obtain ap{u,s). The errors are smaller than for the 
running coupling, so we are able to combine both constant and linear continuum extrapolations 
to control the systematic error from the choice of extrapolation. We estimate the errors using 
a multistage bootstrapping procedure. 

We plot ap in Figure El where the statistical error is in black, and the systematic error 
arising both from the choice of the number of fit parameters and the continuum extrapolation 
is in grey. We also plot the 1-loop perturbative prediction for ap. Our results are close to the 
1-loop prediction, with the running becoming a little faster at strong couplings. 

In the vicinity of an IRFP, we can define an estimator ■j{u) given by 



7(m) 



log I (Tp(m,s) 
log |s| 



:iii 



which is equal to the anomalous dimension 7 at the fixed point 18|], and deviates away from 
the fixed point as the anomalous dimension begins to run. We plot this estimator in Figure JH 
Again the black error bars show the statistical errors, and the grey the systematic errors. We 
see that j{u) is small in most of the range of couplings that we measure. However, it becomes 
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t 


L/a 




8 




10 


10| 


12 


14 


15 


16 


params 
xVdof 


6 
1.79 


5 
0.86 


5 
1.09 


5 
0.62 


5 
0.99 


5 
1.60 


5 
1.63 


5 
1.22 



Table 6: Interpolation best fit parameters for Zp. 



t 


L/a 




8 


9-^ 


10 


lOi 


12 


14 


15 


16 


params 
xVdof 


7 
2.09 


6 
0.46 


6 
1.03 


6 
0.43 


6 
1.18 


6 
0.93 


6 
1.32 


6 
1.47 



Table 7: Interpolation next-best fit parameters for Zp. 

larger at our strongest couplings. Our data is consistent with it reaching values 7 f« 1 that are 
interesting for models of technicolor, although our error bars are large and it is also possible 
that it is as small as 0.135, our lower bound at 'g^ = 4.02, the smallest coupling at which 
a fixed point is consistent with our results using a constant continuum extrapolation. The 
highest value compatible with our data is 7 = 1.03 at 'g^ = 5.52, the highest coupling at which 
we have results for all L. 



5. Conclusions 

In this work we have calculated the running of the Schrodinger Functional coupling ^^ and 
the mass in the continuum limit of SU{2) lattice gauge theory with six flavours of fundamental 
fermions, over a wide range of couplings up to 'g'^ ~ 5.5. 

Our results for the running of the coupling have relatively large errors. This is due to 
the difficulty of measuring the small difference in the coupling between two nearby scales, 
a problem that becomes particularly acute as we approach a possible fixed point where the 
difference falls. We observe that the running of the coupling is slower than the (already slow) 
one-loop perturbative prediction. Our results using a constant continuum extrapolation are 
consistent with the presence of a fixed point above g^ = 4.02, but it is also possible that 
there is no fixed point in the range of couplings we have measured. There is an additional 
uncertainty arising from the choice of continuum extrapolation. 

Our results for the running of the mass are clearer. We find the anomalous dimension is 
small throughout most of the range of couplings we measure, but it becomes larger for our 
strongest couplings, with a possibility that it reaches values around 1. If true, this would be 
very interesting for technicolor models. 

The value of 7 at the fixed point can be predicted using the all-orders beta function 
conjecture 42|. This gives an exact prediction in terms of group-theoretical factors only. If 
the present theory is inside the conformal window, the prediction is 7 = 5/3, which lies outside 
the range measured in this study. Unfortunately, given the uncertainty on the existence of 
the fixed point, any conclusion on the validity of the all-order beta function conjecture is 
speculative at present. 

The accuracy of our results would be improved in particular by using larger lattices, which 
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Figure 4: j{u) using both a constant continuum extrapolation of the two points closest to the continuum, and 
a linear continuum extrapolation. Statistical error using the optimal fit parameters with a linear continuum 
extrapolation in black, systematic error including the choice of continuum extrapolation in grey. 
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would give a larger range of a/L for the continuum extrapolations. This would help to clarify 
the existence and location of the fixed point, and to reduce the errors on the anomalous 
dimension. Calculations to improve the statistics and use larger lattice sizes are ongoing. 
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